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Abstract 
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Magnetization switching in highly anisotropic single-domain ferromagnets 
has been previously shown to be qualitatively described by the droplet theory 
of metastable decay and simulations of two-dimensional kinetic Ising systems 
with periodic boundary conditions. In this article we consider the effects of 
boundary conditions on the switching phenomena. A rich range of behaviors 
is predicted by droplet theory: the specific mechanism by which switching oc- 
curs depends on the structure of the boundary, the particle size, the tempera- 
ture, and the strength of the applied field. The theory predicts the existence 
of a peak in the switching field as a function of system size in both systems 
with periodic boundary conditions and in systems with boundaries. The size 
of the peak is strongly dependent on the boundary effects. It is generally 
reduced by open boundary conditions, and in some cases it disappears if the 
boundaries are too favorable towards nucleation. However, we also demon- 
strate conditions under which the peak remains discernible. This peak arises 
as a purely dynamic effect and is not related to the possible existence of mul- 
tiple domains. We illustrate the predictions of droplet theory by Monte Carlo 
simulations of two-dimensional Ising systems with various system shapes and 
boundary conditions. 

PACS Number(s): 75.50.Tt, 75.40.Mg, 64.60.Qb, 05.50.+q. 
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I. INTRODUCTION 



The next generation of high-density magnetic recording media should have much higher 
storage densities without sacrificing large coercivity. The experimentally observed non- 
monotonic dependence of the coercivity on particle diameter (see, e.g., Ref. 1) implies that 
the optimum size of ferromagnetic grains for use in recording media is small enough to be 
single-domain but large enough to be nonsuperparamagnetic. 2 ' 3 The purpose of this paper 
is to extend earlier studies 4-8 which used Ising systems to model the nonequilibrium sta- 
tistical mechanics of magnetization reversal. Specifically, we investigate the dependence of 
the switching field, which is usually measured in static or slowly increasing fields, on the 
boundary conditions. We show that the boundary conditions strongly influence switching 
phenomena for weak applied fields and small systems, and that this influence is observable 
in the switching field. Some preliminary results of the present study have been published in 



The simplest theory of magnetization reversal in single-domain ferromagnets is due to 
Neel 10 and Brown. 11 In order to avoid an energy barrier due to exchange interactions between 
atomic moments with unlike orientations, Neel-Brown theory assumes uniform rotation of all 
the atomic moments in the system. The remaining barrier is caused by magnetic anisotropy, 2 
which may have contributions from both the local crystalline environment and from the 
macroscopic shape of the particle. Anisotropy makes it energetically favorable for each 
atomic moment to be aligned along one or more "easy" axes. Because the magnetization of 
the entire system rotates uniformly, the free-energy barrier AF separating the metastable 
(magnetization antiparallel to the applied field) phase from the stable (magnetization parallel 
to the applied field) phase is proportional to the volume of the particle. The lifetime r of 
the metastable phase is related to AF by the Van't Hoff-Arrhenius 12 equation 



where f3~ l = k^T is the temperature in units of energy. In what follows, we use units in 
which — 1. Another specific prediction 3 from the uniform rotation model is that the 
switching field increases with system size as 



for a three-dimensional particle. Here H is the asymptotic value of the switching field; 
N is the particle volume; N P is a volume depending on the anisotropy, the waiting time 
r, the saturation magnetic moment of the particle, and the temperature; and Amd is the 
volume above which the equilibrium magnetic structure of the particle consists of two or 
more domains. For A> Amd 5 the switching field decreases, since magnetic reversal can take 
place by means of domain-wall motion. 



Ref. 9. 



A. Micromagnetic theories of magnetization reversal 



r oc exp(/3AF) , 



(1) 




(2) 
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Detailed descriptions of both the static and dynamic properties of fine ferromagnetic 
grains have typically been formulated from micromagnetic models, 13 of which Neel-Brown 
theory is a particularly simple case. This method involves coarse-graining the physical 
lattice onto a computational lattice and then solving the partial differential equations for 
the evolution of the coarse-grained magnetic structures. Although micromagnetics provides 
a good treatment for the anisotropy and demagnetizing fields, 14 it treats thermal effects 
rather crudely, usually just by making the do main- wall energy temperature-dependent. A 
somewhat better approximation for thermal fluctuations within the underlying differential 
equations is to include a Langevin noise term. 15 Yet another possibility is a mean-field type 
approach to switching behavior. 16 A better treatment for thermal and time-dependent effects 
is Monte Carlo simulations (see, e.g., Refs. 4-8, 17-22). Even when the physical phenomena 
can be accurately simulated, however, it will be difficult to understand the results without 
an adequate theoretical basis. It is the purpose of this paper to extend that theoretical basis 
and support it with Monte Carlo simulations for the Ising model with a variety of boundary 
conditions. 



B. Droplet Theory of Magnetization Reversal 

For highly anisotropic systems, magnetization reversal occurs not by means of uniform 
rotation, but by the nucleation and growth of nonequilibrium droplets of the stable magnetic 
phase. 23 The coupling between the magnetic field and the magnetization within a droplet 
favors growth of the droplet, but the interface tension of the droplet favors shrinkage. A 
droplet for which these two forces are balanced is termed a critical droplet. 

Kinetic nearest-neighbor Ising models form a class of highly anisotropic model systems 
that have become popular subjects for Monte Carlo simulations of metastable decay (see 
Ref. 24 and references cited therein). The reasons for this are easy to understand. A 
great variety of exact results have been obtained for the equilibrium two-dimensional Ising 
model in zero field. 25-29 As a nontrivial model for which a number of exact results are 
known, it has become a favorite of researchers in statistical mechanics, and it has been a 
standard in investigations of universality, 30 finite-size scaling, 31 and various approximation 
schemes, such as series expansions. 32 Kinetic Ising models with a variety of stochastic spin- 
flip dynamics have proven useful as simple but nontrivial models for the study and testing 
of ideas in nonequilibrium statistical mechanics, such as dynamic universality and critical 
exponents, 33,34 and metastability. 35 

For these reasons, kinetic nearest-neighbor Ising models have been extensively studied 
as prototypes for switching dynamics in highly anisotropic systems. Square- and cubic- 
lattice Ising systems with periodic boundary conditions have been used to study grain- 
size effects in ferroelectric switching. 36,37 Magnetization reversal in elongated ferromagnetic 
particles has been studied with a one-dimensional model, 38 and a triangular-lattice Ising 
model with mean-field magnetostatic interactions has been shown to reproduce well the 
switching dynamics in Dy/Fe ultrathin films. 17 Square-lattice kinetic Ising models with free 
boundaries were simulated by Melin, 18 and Serena and Garcia 19,20 performed Monte Carlo 
simulations on square-lattice Ising systems with free circular boundaries. In earlier work 4 ' 6-8 
we have considered kinetic Ising models with periodic boundary conditions as models for 
magnetization switching in single-domain ferromagnetic nanoparticles. In the present paper 
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we extend our analytical and numerical study to the effects of nucleation at the particle 
boundary. 



C. Phenomena Affecting Magnetism at Boundaries 

It is well known that the physics of surfaces is far more complicated than one might ex- 
pect based on simple cross-sections of bulk material. 39 For instance, atoms near the system 
boundary may move from their bulk-crystalline positions in order to lower the free energy, 
a process known as reconstruction. Also, the chemical environment of the system boundary 
may be very different from the bulk, since the boundary provides sites for possible chemisorp- 
tion. Furthermore, the reduction in symmetry often changes the magnetic anisotropy at the 
boundary. This is clearly seen in some thin films which become perpendicularly magnetized 
if their thickness is less than some critical value. 40 

The various boundary effects can influence each other. 41 For instance, the ferromag- 
netic metal nickel undergoes reconstruction when oxygen, 42 ' 43 nitrogen, 44 ' 45 carbon, 43,46 ' 47 
sulfur, 43 ' 48 or alkali metals 49 chemisorb on its surface. It has been observed both 
experimentally 50 ' 51 and computationally 52 that clean nickel undergoes a surface structural 
phase transition at the Curie temperature of bulk nickel. Likewise, the electronic changes 
which cause bonding in chemisorption can lead to changes in magnetic properties. 41 For 
example, an oxygen atom adsorbing on Ni(100) reduces the local magnetic moments on the 
nickel atoms in its vicinity. 53 

The breaking of translational symmetry by surfaces can also affect magnetic properties at 
surfaces directly, as it disrupts the spin-wave spectrum and other long-ranged excitations. 39 
This effect can be particularly pronounced in nanoscale particles, since the ratio of "bound- 
ary" atoms to "bulk" atoms can be significant. 54 

As a result of the wide variety of physical phenomena which can affect magnetism at the 
surfaces of single-domain grains, it is clear that an adequate understanding of magnetization 
switching in real ferromagnetic grains requires an understanding of the effects of fairly general 
boundary conditions. In this study we therefore consider in detail the effects of several 
specific modifications to the system boundary. 



D. Experimental Observations of Nanoscale Magnets 

Fine ferromagnetic grains have been studied for many years. For instance, in Ref. 1 
Kneller and Luborsky measured the coercivity and remanence of iron and iron-cobalt par- 
ticles. They found that the coercivity measurements agreed well with Eq. (2) for particles 
sufficiently small to be single-domain, and that the coercivity rapidly drops to the small 
bulk value for multi-domain particles. Until recently, however, such particles could be stud- 
ied experimentally only in powders, which made it difficult to differentiate the statistical 
properties of single-grain switching from effects resulting from distributions in particle sizes, 
compositions, and local environments, or from interactions between grains. 

Recently a variety of techniques have become available which permit one to resolve the 
magnetic properties of isolated, well-characterized, single-domain ferromagnets. 55 One such 
technique is magnetic force microscopy (MFM). 56 Several experiments were important for 
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inspiring our interest in a possible connection between magnetization switching in single- 
domain magnetic grains and metastability in kinetic Ising models. 

In Ref. 57, Chang et al. used MFM to study the switching field of isolated barium ferrite 
fine particles. The shapes of the particles (determined by atomic force microscopy) were 
used to determine the demagnetizing field for each particle, which was subtracted from the 
observed switching field. The resulting effective field has a peak with respect to the particle 
diameter somewhat similar in appearance to that which was reported in Ref. 1, although 
in this case it is believed that the grains remain single-domain, since the peak occurs at a 
grain diameter of about 55 nm. Thus it is difficult to explain the peak in switching field 
reported in Ref. 57 in terms of Neel-Brown theory. In Refs. 4 and 6 we showed that a 
purely dynamic crossover in the switching mechanism can cause such a peak in systems 
with periodic boundary conditions. A central aim of the present study is to understand how 
more realistic boundary conditions affect this peak. 

MFM measurements on individual iron particles 58 showed that as the particle diameters 
increase from 20 nm to 70 nm, the switching fields decrease. The angular dependence of 
the switching field indicates that the magnetization reversal is not coherent at small field 
angles. 

Size effects somewhat similar to those in barium ferrite particles have also been observed 
in nanoscale single-domain Ni bars. 59 The switching field of an isolated bar initially increases 
with increasing bar width (the length of bars used in the experiment was fixed), then reaches 
a maximum at a width of 55 nm, and decreases with further bar-width increase. As the 
authors of that study note, the decrease for wide bars probably occurs because the bars 
change from single-domain to multidomain with increasing width, since only bars of width 
less than 150 nm are single-domain. However, one may speculate that the switching-field 
peak and the decrease in the 50 — 100 nm range could originate from dynamic effects such 
as those found in computer simulations of Ising-like systems. 4 ' 6 

The magnetization reversal mode in long Ni columns was studied in Ref. 60. A weak 
dependence of the switching field on the radius of the column was found, which is not 
consistent with the curling mode. Nucleation at the particle ends was suggested to be 
responsible for this behavior. Measurements of switching fields and their histograms for 
similar Ni wires with smaller diameters ( < 100 nm) give evidence that the magnetization 
reversal in narrow wires results from a nucleation and growth processes. 61 

Peaks in the switching field with respect to the particle diameter have typically been 
attributed to a crossover to a multidomain initial condition. However, if a field significantly 
stronger than the demagnetizing field is initially applied parallel to the particle magnetiza- 
tion, the probability of even a very small domain existing as an initial condition when the 
field is reversed can be made quite small, even for particles that are multidomain in equilib- 
rium at zero field. It is possible that these conditions were achieved in a recent experiment 
by Yang et al., 62 in which they applied a 2T initializing field to nanocrystalline powders of 
the ferrimagnet LiFe 5 8 . Particularly intriguing is the fact that the ratio between the peak 
value of H sw and the value for the largest particle size they report (860 nm) is approximately 
1.7. By comparing the dependence of the lifetime on the exponential part of the nucleation 
rate, it can be shown that for a three-dimensional periodic system that ratio should be 2, at 
least in the limit of long waiting times. 8 It is worth pointing out that surface effects are a 
major concern of Ref. 62, so it may be that surface nucleation is sufficiently discouraged to 
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cause the particles to be similar to periodic systems. Unfortunately, the large polydispersity 
in their powder samples makes detailed interpretation difficult. Microscopic measurements 
of isolated, well-characterized particles under similar conditions would therefore be highly 
desirable. 

Another set of experiments involved the measurement of the probability P that the 
magnetization in 7-Fe203 63 ~ 65 single-domain particles is switched by an opposing field. The 
probability P was measured as a function of applied field for a constant waiting time, and 
as a function of time for a constant applied field. The experiments show that the mean 
switching time (or lifetime) r depends strongly on the strength of the applied field H, which 
indicates that droplet nucleation is occurring — if switching were taking place by domain 
wall motion, r would depend on H much more weakly. 

Magnetization switching in nanoscale particles is often driven by a desire to investigate 
macroscopic quantum tunneling. In a series of papers, 66-69 Wernsdorfer and co-authors 
studied the dynamics of magnetization reversal at very low temperatures, which perhaps 
span the temperature regimes of macroscopic quantum tunneling and thermal activation. 
They measured switching fields and their probability distributions, as well as the probability 
for switching at constant applied magnetic field in cobalt nanoparticles. 68 They suggested 
that the magnetization switching is triggered by a nucleation process and that the reversal 
mechanism in one small particle is more complex than the model of uniform rotation. They 
also observed that the topological fluctuations of the magnetization state, due to sample 
imperfections, affect the switching properties at low temperatures. 67 ' 69 

For some materials, such as thin films of Fe on a sapphire substrate, patterned single- 
crystal islands can be made small enough that they will be single domain, and interactions 
with the substrate can yield a single easy axis. In fact, these islands remain single domain 
over a wide range of sizes and shapes, and their easy axes lie in the same directions indepen- 
dently of the island shape. 70 Consequently, the Ising model is a reasonable first description 
of such materials. 

The experimental evidence discussed above indicates both that nucleation is an impor- 
tant aspect of magnetization switching in nanoscale ferromagnets, and that heterogeneous 
nucleation at boundaries and defects plays at least as important a role as homogeneous 
nucleation in the particle bulk. It provides the motivation for the extension of our previous 
work on bulk nucleation to consider nucleation at the particle boundaries, which we present 
in this paper. 

E. Organization of the Article 

The organization of the rest of this paper is as follows. In Sec. II we describe the boundary 
conditions and dynamics used. In Sec. Ill we briefly review droplet theory for systems 
with periodic boundary conditions. These ideas, with some extensions, are necessary for a 
proper understanding of droplet theory in systems with restricted geometry. In Sec. IV we 
apply the concepts of droplet theory to systems with external boundaries. We particularly 
concentrate on the dependence of the switching field on the particle size and boundary 
conditions. In Sec. V we present our Monte Carlo results and compare them with our 
theoretical predictions. Finally, in Sec. VI we summarize our results. 
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II. DESCRIPTION OF THE MODEL SYSTEMS 

A. Boundary Conditions 

In order to facilitate study of boundary effects such as those mentioned in Subsec. IC, 
as well as the effects of a simply truncated lattice, we have modified the square-lattice Ising 
Hamiltonian (see Subsec. II B below) by allowing the exchange interaction and the applied 
field to be modified at the boundary. This type of modification has yielded a variety of both 
numerical results 71 and exact results 72 for the Ising model in equilibrium — particularly for 
the boundary magnetization and the boundary free energy in the absence of a magnetic field 
applied to the bulk. 25 ' 28 We introduce three types of systems. 

i) Square systems of side L with periodic boundary conditions in only the crystal axis y 
direction are referred to as semiperiodic systems; such systems have the topology of a 
cylinder. Quantities related specifically to these system are denoted by the superscript 
□ in order to distinguish them from their counterparts for other types of systems. 

ii) Systems consisting of all sites lying within a circle of diameter L centered midway 
between sites [i.e., at (|,|) if there is a site at both (0,0) and (1,1) on a square lattice] 
and with open boundary conditions are referred to as circular systems. We use the 
superscript to denote quantities specific to the circular lattices. 

iii) Systems consisting of all sites lying within an octagon of width L, centered midway 
between sites and with periodic boundary conditions in no direction are referred to 
as octagonal systems. In an attempt to make the boundary more uniform, half the 
nearest-neighbor bonds on the edges of the system parallel to the x- and y axes are 
omitted. See Fig. 1. 



B. The Hamiltonian 

The Ising model is defined by the Hamiltonian 

Ho = -J^SiSj - H^Si , (3) 

where Sj = ±l is the z component of the dimensionless magnetization of the atom (spin) at 
site i, J > is the ferromagnetic exchange interaction, and H is the applied magnetic field 
times the single-spin magnetic moment. The sums an d J2i run over an nearest-neighbor 
pairs and all sites on a lattice, respectively. The dimensionless system magnetization is given 
by 

m = N- 1 J2s i , (4) 

i 

where N is the total number of sites in the system. The lattice constant is set to unity. In 
this article we consider different shapes and boundary conditions for the system, but the 
underlying lattice is always taken to be square. For the reader's convenience, we note that 
the critical temperature of the two-dimensional Ising model is 2 J/ log^ + v^) = 2. 26919 J. 25 
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The selection of the Ising model is equivalent to requiring a very large (infinite, in fact) 
anisotropy constant. Although magnetic materials used in magnetic recording media require 
comparatively large anisotropy constants, 73 the microscopic anisotropy tends to be much 
smaller than the exchange energy. However, in some applications, such as many thin films, it 
is convenient to use Ising spins to represent individual grains that are superferromagnetically 
coupled to make up the system (see, e.g., Refs. 17, 74, 75). If these coupled grains reverse 
their magnetization through uniform rotation, as in Neel-Brown theory, 10,11 the anisotropy 
barrier for a grain is the product of the anisotropy barrier for a single atom and the grain 
volume. Thus, although the present work is intended as a step towards a quantitative 
microscopic theory for nanoscale ferromagnetic particles, it can also be used to describe 
systems consisting of superferromagnetically coupled grains. 

In order to give greater flexibility in dealing with boundaries, we modify the Hamilto- 
nian in Eq. (3) by allowing an additional field i?s only at the boundary and an additional 
coupling J 2 that connects only nearest-neighbor boundary spins (see Fig. 1). This yields 
the Hamiltonian 

= T-Co ~ s i s i - HzJ2 s i ■ ( 5 ) 

The sums an d Si E run over all nearest-neighbor pairs of boundary sites and over all 

boundary sites, respectively. 



C. Dynamics and Simulation Methods 

The relaxation kinetics are simulated by the single-spin-flip Metropolis dynamic 76 with 
updates at randomly chosen sites. The acceptance probability in the Metropolis dynamic 
for a proposed flip of the spin at site a from s a to — s a is defined as 

Wu(sa^-Sa) = min[l, exp(-/3A£ Q )] , (6) 

where AE a is the energy change due to the flip. A rigorous derivation from microscopic 
quantum Hamiltonians of the stochastic Glauber dynamic 77 used in Monte Carlo simulations 
of Ising models has been established in the thermodynamic limit under certain restrictions. 78 
Both the Glauber and Metropolis algorithms are spatially local dynamics with nonconserved 
order parameter (the dynamic universality class of Model A in the classification scheme of 
Hohenberg and Halperin 33 ) and are therefore expected to differ only in nonuniversal features. 

We implement the Metropolis dynamic both by the original, straightforward Metropolis 
algorithm 76 and by the "refusal-free" n-fold way algorithm. 79 ' 80 The n-fold way algorithm is 
much more efficient than the original Metropolis algorithm at low temperatures, where the 
latter requires many attempts before a change is made. 

The switching process is simulated by starting from an initial state fully magnetized 
opposite to the applied field, i.e. all spins are in the state Sj = +1 and the field H < 0. The 
definition of the lifetime of the metastable state is the mean first-passage time to the cutoff 
magnetization m = 0: 

r=(t(m = 0)). (7) 
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From the lifetimes measured for various strengths of the magnetic field and various system 
sizes, we numerically determine the switching field, H sw , which is defined as the absolute 
value of the field required for a system of a given size to exhibit a given lifetime. Thus, the 
switching field is a function of the temperature, the system size, and of the waiting time. 
Note that the terms "lifetime" and "waiting time" both stand for the same quantity, r. We 
use "lifetime" whenever r is treated as a dependent variable t(L, H, T), and we use "waiting 
time" when r is considered an independent variable. 



III. DROPLET THEORY FOR HOMOGENEOUS NUCLEATION 

In this section we briefly review the droplet theory of nucleation for systems with periodic 
boundary conditions. Most of this material is treated in greater detail in earlier works (see, 
e.g., Refs. 4, 6-8, 24, 81, 82) but the concepts are needed when we deal with systems with 
distinct boundaries in the next section, and some of the details we present here have not 
been published before. 



A. General Considerations 

The central problem in nucleation theory is to evaluate the free-energy barrier for nucle- 
ation of the equilibrium phase, AF, which is the free-energy difference between the system 
containing a single critical fluctuation and the same system in the homogeneous metastable 
phase. Once this is determined, the dominant field and temperature dependence of the 
nucleation rate per unit volume is obtained from the Van't Hoff-Arrhenius relation, Eq. (1). 
However, in contrast to the mean-field Neel-Brown theory, AF is not proportional to the 
system volume. 

For the purposes of this paper, it suffices to adopt an approximation in which the system 
is divided into two distinct regions: a metastable background in which the magnetization 
approximately equals the zero-field spontaneous magnetization m sp , and a region in which 
the magnetization is approximately — m sp , parallel to the applied field H <0. This amounts 
to setting the magnetic susceptibility x=0 and corresponds to ignoring Gibbs-Thomson type 
corrections. 83 For details on the validity of this approximation for kinetic Ising systems in 
moderately strong fields, see Refs. 7 and 84. 

First we consider the case that the region of stable phase is a single droplet of radius 
R. For the time being, we do not take into account the fact that the droplet can nucleate 
anywhere in the system. Such entropy-related contributions to the free energy of a droplet 
will be included later. The free-energy barrier corresponding to creating a droplet of radius 
R is then 

F dr (R) = ft [daR d - 1 - 2\H\m sp R d ] , (8) 

were a is the interface tension in a symmetry direction of the lattice and d is the spatial 
dimension. The quantity ft is a temperature-dependent factor which gives the volume of 
a droplet of radius R as VlR d ] thus for two-dimensional systems at sufficiently high tem- 
peratures that o is not very anisotropic, ft ps it. Here and elsewhere in this paper, we set 
the free energy equal to zero in the uniform metastable phase. The maximum of F dr (R) 
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corresponds to a droplet for which the tendency to grow, due to the coupling of the field and 
the magnetization, is balanced by the tendency of the droplet to shrink, due to the interface 
tension. The radius of this critical droplet is 



R c — 



(d-l)a 



(9) 



2\H\m. 



The free-energy cost of the critical droplet is obtained by using R c in F dr (R): 



d-l 



) 



d-l 



AF SD = Va' 



A 



(10) 



2\H\m, 



sp 



The subscript SD stands for single-droplet, as will be explained in the next subsection. 
The free-energy barrier AFsd is a function of the temperature T through Q, a, and m sp . 
However, it is not dependent on the system size. Later it will become clear that this is an 
important point. 

Inserting AF SD into Eq. (1) one gets the dominant contribution to the nucleation rate per 
unit system volume, I(T,H). A field-theoretical saddle-point calculation gives corrections 
due to fluctuations about the critical droplet configuration, 24 ' 85 ' 86 yielding the final result: 



where B(T) is a nonuniversal prefactor. The prefactor exponent K is 3 for the two- 
dimensional Ising model 82 ' 85 ' 86 and —1/3 for the three-dimensional Ising model, 86 assuming 
local diffusional dynamics. 



Scaling arguments discussed in detail in Refs. 24, 81, and 82 reveal four distinct regions in 
the space of field and system size, in which the switching proceeds in qualitatively different 
ways depending on the number of critical droplets which form during the switching process. 
Below we consider separately the three regimes that are relevant to the range of relatively 
weak fields studied in the present paper. 

The coexistence (CE) regime 81 ' 82 is the part of the L-H plane characterized by systems 
which are too small to accommodate a critical droplet. Instead, the free-energy barrier 
separating the stable and metastable phases comes from a "slab" of stable phase, 87 separated 
from the metastable background by two interfaces of area L d ~ x . Since the slab volume can 
be increased simply by separating these interfaces without increasing the interface area, once 
such a slab has formed, it almost always continues to grow. The free energy of a "critical" 
slab is found by equating the interface part of the droplet free energy [the first term in 
Eq. (8)] with the interfacial free energy of a slab, 2aL d ~ 1 , and inserting the resulting value 
of R into Eq. (8). The result is: 



I(T, H) = B(T)\H\ K exp {-/? [AF SD + O (h^)] } , 



(11) 



B. Modes of Metastable Decay 




(12) 
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Below we show that the CE regime corresponds to fields weaker than a crossover field 
proportional to L^ 1 . As a result, AFqe depends on L as L d_1 . The dominant contribution 
to the lifetime for slabs is obtained by inserting AFqe in Eq. (1). In addition there may 
be prefactors analogous to those in Eq. (11), but their forms are not known. We therefore 
obtain the following expression for the lifetime in the CE regime: 

t ce (L, H, T) » A(T, L) exp [(3AF CE \ , (13) 

where the prefactor A(T, L) is nonuniversal. The dominant behavior of this approximate 
result agrees with other studies. 88-90 The switching field H sw was explicitly defined in Sub- 
sec. II C as the absolute value of the field that corresponds to a specified waiting time r for 
given L and T. An approximate form for H sw , valid in the CE regime, is obtained by solving 
Eq. (13) with AF CE given by Eq. (12) for \H\ at fixed r, while ignoring the pre-exponential 
L dependence: 

[ for L < d^L ThSp (r) 

H sw (L,r,T)^\ \Qd d /2 d ] lLj ( \n(r/A)\ i ■ ( 14 ) 

The second line in this equation represents a rapidly increasing function of L. Here 



-^ThSp(T) = 



d Hr/A) 



2[3a 



i 



(15) 



is a value of L which can be shown to correspond to the system size at which the switching 
process crosses over into the single- droplet regime (see below), known as the "thermodynamic 
spinodal" (ThSp). 81 ' 82 The thermodynamic spinodal is the point at which the free energy 
of the critical droplet discussed in Subsec. Ill A above equals that of a slab. It corresponds 
to a first-order phase transition in the fixed-magnetization ensemble. 87 ' 91 Since the free 
energies of both these excitations are proportional to their interface areas, this corresponds 
to dQ,R d ~ 1 2L d_1 . The value of the switching field at the thermodynamic spinodal is 
therefore given by 

1 (ttd\^ (d-l)a 

#Ths P ~^ hH n ( \ m ' 16 

2 V 2 / ( L ThSp(l") - ^Hsp 

where i is a weakly temperature-dependent phenomenological parameter of order unity (see 
Table I), which is fitted to simulation data in order to compensate for finite-size effects. 8 In 
the remainder of this paper we do not consider corrections of the type accounted for by £, but 
we point out that without such corrections, the peak in the switching field found numerically 
for the small systems we can simulate is not in very good agreement with Eq. (16). 

The single- droplet (SD) regime 81 ' 82 is the part of the L-H plane characterized by systems 
for which the first critical droplet to nucleate grows to the size of the system before a second 
critical droplet nucleates. The lifetime in the SD regime approximately equals the average 
time until the first droplet nucleates, 



TSD 



t nuc = (L d l) 1 oc L~ d exp [/3AF SD ] . (17) 
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Solving this equation with AFsb given by Eq. (10) for \H\ at fixed r, we obtain an approx- 
imate form for the switching field in the SD regime: 



H sw (L,t,T) 



2m sp 



fitter 



In (rL d /B) 



for L > L ThSp (r) . (18) 



This is a monotonically decreasing function of L. Due to the neglect of prefactors in the 
derivations of Eqs. (14) and (18), these estimates for the switching field do not match 
perfectly at LihSp- One function of the adjustable parameter £ in Eq. (16) is to provide this 
matching. 

Equations (14) and (18) show that the derivative of H sw with respect to L at fixed r is 
positive in the CE regime and negative in the SD regime. The crossover between these two 
regimes must therefore be associated with a maximum in the switching field. We emphasize 
that the presence of this peak is due solely to a nonequilibrium crossover and not to the 
effects of magnetostatic dipole-dipole interactions 7 which can cause the equilibrium system 
to be composed of more than a single domain. This is in contrast to Eq. (2), which indicates 
that the switching field increases monotonically as long as the system is single-domain in 
equilibrium. 

In both the CE regime and the SD regime, switching is abrupt, with a negligible amount 
of time spent in configurations with magnetizations significantly different from =tm sp , and 
switching is a Poisson process. This phenomenon, in which the entire system behaves as 
though it were a single magnetic moment, is known as superparamagnetism. 2 ' 3 As a conse- 
quence of the Poisson nature of the switching process, the standard deviation of the lifetime 
for an individual grain is approximately equal to the mean lifetime, r. Because of the random 
nature of switching of these two regions of the L-H plane, they have been jointly called 81,82 
the "stochastic" region. 

We can combine Eqs. (10), (12) and (16) (the latter with £=0) into the convenient form 



AF 



AF< 



SD 




< 1 

(19) 

> 1 

ThSp 



which is shown for d=2 in Fig. 2. By comparing Eqs. (9) and (16), we find that the argument 
in this equation has a transparent interpretation in terms of the relative sizes of the system 
and the critical droplet: 



-f^ThSp \£ldj 2R C 

Although Eq. (19) is a monotonic function, it provides a clue to the existence of a peak in 
the switching field at the thermodynamic spinodal. For \H\ >H ThSp (i.e., in the SD regime), 
the free-energy barrier for the nucleation of a droplet at a specified site becomes independent 
of the system size. This is basically because the shape of the droplet is independent of its 
radius. However, the system size still enters into the lifetime through the number of available 
nucleation sites, as shown in Eqs. (13) and (17). These entropic corrections to Eq. (19) and 
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their analogues for systems with nonperiodic boundary conditions are the main cause of the 
peak observed in the switching field. They are discussed in greater detail in Subsec. IV C 
below. 

For fields that are not too strong, the radial growth velocity of a supercritical droplet 
is linear in \H\: v = z/|if |. 4 ' 84,92 ~ 95 This leads to a competition between nucleation and 
growth. The multidroplet (MD) regime 81 ' 82 is the part of the L-H plane where the finite 
growth velocity allows other critical droplets to nucleate before the first critical droplet has 
grown to the size of the system. The crossover between the SD and MD regimes is known 
as the "dynamic spinodal" (DSp). 81 ' 82 A reasonable criterion to locate the DSp is that the 

nucleation time, t nuc = {L d l) , and the time it takes a droplet to grow to a size comparable 
to L should be equal: 4 ' 91 



L d I(T,H) = (2n)- 1/d L- 1 u\H\ 
This yields the asymptotic relation 4 ' 82 



(21) 



DSp 



(d-1) 



2m 



sp 



/3tta a 



(d+l)\nL 



(22) 



However, relatively large systems (L > 10 3 10 4 for d—2) are required for the contribution 
described by Eq. (22) to be larger than the correction terms in the numerical solution of 
Eq. (21) (see Fig. 11 of Ref. 91). 

In the MD regime the magnetization is a self-averaging quantity and develops determin- 
istically with time according to "Avrami's i aw: " 24 > 82 > 96 ~ 98 



m(t) ~ 2m sp exp 



t \ d+1 
tmd/ 



m. 



Sp 5 



(23) 



where 



(d+l) In 2 



i 

" d+l 



oc exp 



SD 



d+l 



(24) 



gives the lifetime. The factor d+l which occurs in Eqs. (23) and (24) is the space-time 
dimension of the system. Note that tmd is independent of L. Details of switching in the MD 
regime for systems with periodic boundary conditions are investigated elsewhere. 84 ' 99,100 

In the figures that show Monte Carlo results for the switching field versus L for various 
boundary conditions (Figs. 7, 9, 11, 12), the corresponding values for periodic boundary 
conditions are shown as dotted curves. The general shape expected from the above discussion 
— zero for small L, then a sharp rise, followed by a decline in the SD regime towards a plateau 
in the MD regime — is clearly observable. 



IV. DROPLET THEORY FOR HETEROGENEOUS NUCLEATION 

In this section we examine the free-energy barriers between the metastable and stable 
phases for two simple system geometries described in Subsec. II A: semiperiodic and circular 
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systems. In doing so, we restrict our theoretical considerations to two-dimensional systems, 
as we also do for our numerical study. However, the results can be generalized to d=3 in 
a straightforward manner. For simplicity, we also ignore the anisotropy of the interface 
tension. This is an excellent approximation for intermediate and high temperatures. In 
what follows, whenever we refer explicitly to quantities that correspond to homogeneous 
nucleation, such as cross-over fields or free-energy barriers, the special case d—2 and Q=tt is 
to be understood. We also maintain a nomenclature such that the terms "volume," "area," 
and "length" correspond to the two-dimensional equivalents of (length)^, (length)^" 1 , and 
(length) 1 , respectively. 

As with systems with periodic boundary conditions, the dominant contribution to the 
lifetime is an exponential dependence on the free-energy barrier for a critical fluctuation cen- 
tered at a particular site, as shown in Eqs. (13), (17), and (24). In Subsecs. IVA and IVB 
we obtain the corresponding free-energy barriers for semiperiodic and circular systems, re- 
spectively. 

As was also discussed in Sec. Ill, the dominant behavior of the lifetime is modified by 
corrections involving the number of available nucleation sites and the applied field. We show 
in Subsec. IV C that these corrections can be viewed as entropy corrections to the free-energy 
barrier. The existence, size, and location of the peak in the switching field are determined 
by the manner in which these corrections depend on H and L. 

The surface of a droplet of stable phase forming on the system boundary consists of two 
distinct parts. The first is an interior interface separating the droplet from the metastable 
background. This part of the droplet surface contributes a term to the total free energy of 
the system which is the product of its area and the interface tension, a. The second is the 
part of the system boundary which is also part of the droplet. This exterior droplet surface 
contributes a term to the total free energy of the system which is the product of its area and 
the "boundary tension," <j£. This boundary tension expresses the difference in free energy 
per unit area between the system boundaries limiting the two different bulk phases. 

For a two-dimensional system, we can exploit the exactly known free-energy density of 
the one-dimensional Ising model 101 ' 102 in an external field h, 

fi(h) = -T In je^- 7 cosh[/3/i] + y/e~^ J + e 2 ? J sinh 2 [/3/i]| , (25) 

to obtain a rough approximation for the boundary tension. Consider the chain of boundary 
spins. In addition to the external field H, they are also subject to interactions with the 
interior of the system, which we take into account in a mean-field approximation. This 
results in an additional effective field ±Jm sp , the sign of which depends on whether (+) or 
not (— ) there is a droplet attached to the boundary of the system. In this way, one can get 
the contribution to the free energy from the boundary layer of spins as the free energy of 
a one-dimensional Ising system in an appropriate field. Then, an estimate of the boundary 
tension is given by 

<r s » [/:(# + + Jm sp ) - h{H + H s - Jm sp )] - [f^H + Jm sp ) - fj(H - Jm sp )] . (26) 

Here, the first and second terms represent the free-energy difference from the metastable 
phase due to the external boundary of the droplet with and without the applied surface 
field H%, respectively. Note that er s = when if s = 0. This estimate is less accurate at 
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temperatures closer to the critical point, where the mean-field nature of the approximation 
causes a-£ to be overestimated. 

Having estimated the boundary tension, we obtain the contact angle 9 Y of droplets 
forming at the system boundary from Young's relation: 103 



For <t 2 > a, Eq. (27) yields 9y = ir, corresponding to a completely "dry" boundary For 
°"s — 0, one has 9y = vr/2, and for < —cr, 9y = so that the stable phase completely 
"wets" the boundary. We emphasize that although as depends on H, this dependence is 
practically negligible in the region of interest for this study. This is documented in Fig. 3, 
where we have plotted 9y as a function of H% for two different values of H . 

The determination of the droplet shape is directly related to studies of wetting. 104 Some 
readers may be aware that a wetting problem has been solved exactly by Abraham 105 for 
the two-dimensional Ising model, and that the droplet shape in that case is an ellipse, 106 
not a circle. However, the systems Abraham studied were in zero field, and so this result 
and other zero-field studies 107 are not directly applicable to the current topic. 



As is the case for periodic boundary conditions, the decay rate is controlled by the 
nucleation process that has the lowest free-energy barrier. In addition to the bulk droplet 
and the slab connecting two sides with periodic boundary conditions, which were considered 
in Sec. Ill, we now also must consider a droplet nucleating on one of the open system 
boundaries [Fig. 4(a)] and a slab configuration formed by an interface parallel to the open 
boundaries [Fig. 4(b)]. Both these configurations in general give lower free-energy barriers 
than the ones previously considered. We therefore assume for the moment that nucleation 
on the boundary is preferred over nucleation in the bulk, and we neglect entropic terms 
which arise from translational invariance and other corrections. These issues are addressed 
in Subsec. IV C. 

The free energy of a system composed of two "slabs" (si) with magnetizations ±m sp 
[Fig. 4(b)], corresponding to a system magnetization m, is 



The superscript □ is used to designate quantities specifically related to the nonperiodic 
boundaries of a semiperiodic system. 

The droplets form preferentially on an open boundary with a contact angle 9y given by 
the Young equation, Eq. (27) [see Fig. 4(a)]. The magnetization of a system containing a 
single droplet of radius R is then given by 



ctcos(#y) + o"s = . 



(27) 



A. Semiperiodic Systems 



F°(m) = L(a + <j s ) + L 2 \H\ (m - m sp ) . 



(28) 




(29) 



where 



f(0v) 



29 Y - sin(2fl Y ) 
2^ 



(30) 
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The droplet covers a boundary area equal to 2ffsin(#Y) and has an interior interface of area 
2R9 Y . The difference from the free energy of a uniform metastable system is given by 



F d D (m) = 2Rsm(6 Y )ax + 2R9 Y a + L 2 \H\ (m - m sp ) 



= 2tt/(0 y ) (R<j — R 2 \H\m 



sp 



= aLJ2irf(0 Y ). l - — + L 2 \H\ (m - m sp ) . 



(31) 
(32) 



Equating Eqs. (28) and (32) yields the "droplet-slab" (ds) magnetization, at which the 
droplet and slab free energies are equal, 



m 



ds 



m. 



sp 



1 



(1 — COS^y) 

2nf(9 Y ) 



(33) 



Using the condition (d/dm)F^(m)\ m n = 0, we find the thermodynamic spinodal (to leading 



order in 1/f ), 



-"ThSp 



2/(0 



Y 



7TO" 



2/(0 



Y 



cos^y 2Lm, 



sp 



COS0^ 



-if 



ThSp 



(34) 



Note that if ThSp < -^ThSp for all 0y < ir. 

To evaluate the free-energy barrier for \H\ > H^ hS , we find the critical-droplet radius 
from Eq. (31) by setting (d/dR)F^ = 0. This yields Eq. (9), just as for periodic boundary 
conditions. Thus the critical droplet radius is independent both of L and of Y . The free- 
energy barrier is given by F^(R C ): 



Aif D 



f(0y) 



1T/J 



2 if m. 



/(0 Y )AF SD 



(35) 



sp 



Note that AF s n D < AF SD for all 9 Y < n. For |if| < ffxhs P > the free-energy barrier is just 

/•;tV4): 



AF C D E 



La(l — cos0 



Yj 



L 2 \H\m. 



cos0 Y ) 2 



sp 



2tt/(0 y ) 



(36) 



It is worth pointing out that the above results for the special case of 9 Y = ir (i.e., a completely 
"dry" boundary) are identical with the results for periodic boundary conditions with d = 2 
and Q = n. In particular, Eq. (31) corresponds to Eq. (8), Eq. (34) corresponds to Eq. (16), 
Eq. (35) corresponds to Eq. (10), and Eq. (36) corresponds to Eq. (12). 

Just as for periodic boundary conditions, we can combine Eqs. (34), (35) and (36) to 
obtain 



AF n 



\H\ 



\H\ 



ttU 

.-"ThSp, 



ttU 

,-"ThS P/ 



for 
for 



I if I 



TTU 

-"ThSp 
Iffl 



< 1 



(37) 



TTU — 

-"ThSp 



> 1 
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which is the equivalent of Eq. (19) for d=2 and with the 8y dependent H^ hSp instead of 
-^ThSp- This relation is illustrated in Fig. 2. As in Eq. (19), the argument has a simple 
interpretation in terms of the relative sizes of the system and the critical droplet: 

\H\ l-cos6 Y L (3g) 



B. Circular Systems 

Consider a circular droplet growing at the boundary of a circular system with a wetting 
angle of 6, illustrated as 6y in Fig. 5. If the droplet has radius R, then it subtends an arc of 
the particle boundary corresponding to an angle 20, where is given by 

(2^0) = ArcTan^A +C os(0),sin(0)) . (39) 



The two-argument function ArcTan(x, y) is the angle defined by 



108 



X 

cos (ArcTan(x, y)) = . = ,. , (40a) 

Vx 2 + y 2 



sin (ArcTan(o;, y)) = y , (40b) 

y>X 2 + y 2 

and 

< ArcTan(x, y) < 2vr . (40c) 

To simplify the notation, we omit the arguments of the function in what follows. The 
volume of the droplet can be expressed as 

V= ^ + {e -^-^p?l, (41) 

and the total free energy is given by 

F(R, 9) = Lcjxj^ + 2Ra(9 - 0) - 2\H\m sp V . (42) 

In order to obtain the wetting angle 6 we minimize T while holding V constant, so that 

dT_dV_ _ dTdV_ 

After factorization, this condition reads 



. (43) 



as + a cos(#) (2R + L cos(0) ) ( - 9 ) + L sm{9) 



ALR 2 



L 2 + 4R 2 + 4RL cos(fl) 



= . (44) 
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It is clear that among its solutions, only the solution of the first factor in the numerator, 
corresponding to Young's relation, i.e., 9 = 9 Y , is acceptable. This is because that is the 
only solution which depends on a. Using this we find the free energy to be given by 

F°(R)=F(R,9 Y ) 

= -L(f)a cos(# Y ) + 2Ra(9 Y - 0) - 2\H\m sp V . (45) 

The superscript is used to designate quantities specifically related to the boundary of a 
circular system. 

Note that, in contrast to systems with periodic boundary conditions in at least one 
direction, there is no slab solution for circular systems. However, the droplet free energy, 
F e (R), depends quite differently on R and L in the two limits, L/2R— >0 and L/2R— >oo. 
In the former case, which corresponds to a droplet with an almost flat inner interface, it is 
easy to see from Eq. (39) that — > 9 Y . In this limit, F e (R) is therefore independent of R 
and linear in L. In the latter case, which corresponds to a droplet which is so small that the 
system boundary appears flat in comparison, we expand in powers of 2R/L: 



v hi / 9 A? \ 

~ Jl sin (0) _ \ t±\ cc . s (fl) S in(0) + . . . . (46) 

It is easy to see that in the limit L/2R-^oo, Eq. (45) with V given by Eq. (41) is identical to 
Eq. (31) for semiperiodic systems, as expected. Thus, in this limit the free energy is linear 
in R and independent of L. 

A simple calculation shows that the critical radius R c is not affected by the boundary 
conditions and is again given by Eq. (9). As a result, we obtain the free-energy barrier 

AF G = F & (R c ) . (47) 

Since lim£/2.R c ->oo AF Q = AF§ D , it is reasonable to try to express the crossover of AF & 
with L/2R C in a scaling form similar to Eqs. (19) and (37). Specifically, we get 

AF 9 Y -^x,9 Y )[l+x 2 + 2xcos(9 Y )]+xsm(9 Y ) 

, (48a) 



AF s n D 7rf(9 



Yj 



where 



x- L - n _l-cos# Y \H\ ^ gb ^ 

2R C 2 H ThSp tt/(#y) -f^ThSp 

and i^ThSp ^ s gi ven by Eq. (34). Note that x does not depend on 9 Y ; it simply expresses the 
system size in units of the critical droplet size. Equation (48) thus shows that the free-energy 
barrier for circular systems, when normalized by the free-energy barrier for a semiperiodic 
system with the same value of 9 Y) is a function of the field only through \H\/H^ hSp . How- 
ever, the absence of a slab solution implies that there is no sharply defined thermodynamic 
spinodal for circular systems. The crossover is therefore much more gradual than in sys- 
tems with periodic boundaries in at least one direction. Figure 2 shows AF & / AF$ D for two 
different values of 9 Y . 
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C. Entropy Effects Due to Translational Invariance 



As mentioned briefly in Subsecs. Ill B, IV A, and IV B above, the scaling relations for the 
free-energy barriers [Eqs. (19), (37), and (48) for periodic, semiperiodic and circular systems, 
respectively], if interpreted as scaling relations for the logarithm of the metastable lifetimes, 
do not predict a peak in the switching field when considered as functions of the system size 
at constant waiting time. This peak, the dependence of which on the boundary conditions 
is the main topic of the present paper, is chiefly the result of entropic terms that correspond 
to the number of available nucleation sites. When these terms are considered as corrections 
to the free-energy barriers of the droplets discussed above, they lead to corrections that are 
logarithmic in L. They decrease the free-energy barrier as the system size increases, and, as 
a consequence, the switching field is a decreasing function of the system size. We note that 
these entropy effects are analogous to the translational corrections to the nucleation rates in 
fluids, previously introduced by Lothe and Pound. 109 This argument applies only in the SD 
regime, and the switching-field maximum can be observed only if the entropic terms start 
to be discernible before the crossover into the MD regime sets in. 

Other terms in the free-energy barrier that are logarithmic in the applied field arise from 
power-law prefactors in the nucleation rate, such as the one explicitly included in Eq. (11). 
These terms are identical for droplets nucleating in the bulk and droplets nucleating on the 
system boundary, at least within the approximations described in this paper. The growth 
times also need to be considered. These yield terms proportional to \n(L/\H\), as well 
as constant terms. Only the constant terms are different for bulk and surface droplets. 
Consequently, none of these terms are important in the lowest-order approximation. We 
shall find it more convenient to consider the terms arising from power-law prefactors in the 
nucleation rates part of the free-energy barriers, but to treat the growth time separately. 

We begin with a system with semiperiodic boundary conditions for which switching 
occurs through nucleation at the system boundary, which is of area 2L. Then Eqs. (11), 
(17), and (35) yield for the lifetime: 110 

2 

TlnCc = AFsd = f(&y)^m Tln ( 2L ) " TKln \ H \ ■ ( 49 ) 

Z\H \m sp 

Here and henceforth, the notation AF is taken to mean a free-energy barrier for a droplet 
in which the entropy-related corrections (or, in other words, the "effect of the number of 
nucleation sites") as well as terms proportional to Kin \ H\ are included. Similarly, the free- 
energy barrier for nucleating a single droplet in the bulk is obtained from Eqs. (10), (11), 
and (17): 

TlnW = AF SD = ^ Tln(L 2 ) - TK\n \H\ . (50) 

Z\n\m sp 

Nucleation in the bulk will be^the dominant decay mode only if there exists a range of 
L and H such that AF^d < A.F SD . The crossover point in the L-H plane can therefore 
be found as a function of the nucleation time by solving Eqs. (49) and (50) together. For 
sufficiently weak fields we can ignore the ln|if | terms. If we further make the assumption that 
the waiting time r approximately equals the nucleation time, which is valid for single-droplet 
decay, we obtain: 
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/W [2/(^-1] 
xl ; 2m sp ln(4r) 1 ; 

and 

L x (r)» (2 T [ 1 - / ^) 57 ^ T . (51b) 
These solutions are physically acceptable only if 

5 < /(0y) < 1 , (52) 

which means that there can be a transition from boundary nucleation to bulk nucleation 
only if 

| < 0y < vr . (53) 

Eq. (51b) gives an estimate of the system size above which bulk nucleation dominates over 
nucleation at the boundary for a given waiting time. On the other hand, in systems smaller 
than L x (r), the majority of the critical droplets will nucleate at the boundary. We emphasize 
that the crossover between these two regimes is not sharp, but we have observed both regimes 
in Monte Carlo simulations, with a crossover in semiquantitative agreement with Eq. (51). 

The above arguments indicate that for 0y > 7r/2, there exist circumstances in which 
the first droplet to nucleate forms in the bulk, rather than on the boundary. However, one 
may still ask whether the transition from boundary-dominated nucleation to bulk-dominated 
nucleation takes place in the SD regime or in the MD regime. In analogy with Eq. (21), the 
dynamic spinodal for semiperiodic systems may be defined by the equality of the nucleation 
time and the growth time: 

L exp (-/3AFg D ) = it- l/2 L- 1 v\H\ , (54) 

with the asymptotic solution 

^-^sSo-i^-*- (55) 

In order for bulk nucleation to dominate surface nucleation in the SD regime, we must have 
i^Dgp > H x , which yields 

f(0y)>l- (56) 

We designate as #tmik the value of 9y which satisfies Eq. (56) as an equality, and we find 
numerically that #buik~0.5857r. 

In summary, the changes in the decay mode from boundary-dominated single-droplet 
decay with increasing L proceed differently, depending on the value of 9 Y . 

• For #y~71"/2, nucleation at the boundary is always more likely than nucleation in the 
bulk. As L increases, the system crosses the boundary dynamic spinodal, H^ Sp , and 
the decay changes from the boundary-dominated single-droplet mode to the boundary- 
dominated multidroplet mode. Bulk nucleation does not occur for any value of L. 
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• For 7r/2<#Y~#buik, multidroplet nucleation at the boundary becomes dominant at a 
value of L<L X . As L is increased further, bulk multidroplet nucleation becomes the 
dominant decay mode. 

• For #buik~#Y~7r, bulk single-droplet nucleation dominates for L>L X . For sufficiently 
large L the system crosses the bulk dynamic spinodal, Hvs p , and the bulk multidroplet 
decay mode becomes dominant. 

• For 6 , y = 7t, bulk nucleation is dominant for all L. 

The effects of the entropic terms on the switching field are shown in Fig. 6. 

Although here we only explicitly consider the effects of translational invariance for 
semiperiodic systems, at high temperatures we expect no significant difference for systems 
with fully open boundaries, such as circles. At low temperatures, the anisotropy of the in- 
terface tension may cause some nucleation sites to be preferred on the boundaries of circles, 
thus diminishing the translational invariance. 

V. MONTE CARLO SIMULATIONS 

The theoretical results of the previous section deal most carefully with those effects which 
should be dominant in the limit of long waiting time, where the switching field is small and 
the droplets are large. In this section we use Monte Carlo simulations to demonstrate that 
the results of the previous section hold quite well even when the lifetime is not extremely 
long, so that the switching field is rather large and the droplets rather small. 

A. Semiperiodic Systems 

In Fig. 7 we present the switching fields for semiperiodic systems as functions of the 
system size for T = 1.3 J ~ 0.57T C and for different values of the surface field H%. Panels 
(a) and (b) of the figure correspond to the waiting times r = 1000 and 30000 Monte Carlo 
Steps per Spin (MCSS), respectively. In both cases we see similar behavior, depending on 
the value of For large ifs, a semiperiodic system behaves very much like a periodic one. 
The peak of the switching field is pronounced because there is a well distinguished region 
of system sizes in which single-droplet nucleation is dominant. For = 1.0 the critical 
droplet typically nucleates in the bulk even in small systems, whereas for H s = 0.5 it usually 
appears at the boundary. As the system size increases, nucleation in the bulk becomes 
the more important process in both cases, and the switching fields gradually approach 
those measured for the periodic system. For H^, = 0.5, this corresponds to the crossover 
from single-droplet nucleation at the boundary to single-droplet nucleation in the bulk, 
as described for # bulk < 9 Y ^ n in Subsection IV C. This switching-field behavior also 
qualitatively corresponds to the expectation based on our estimates of the contact angles of 
the droplets. From Fig. 3 it is seen that for = 1.0 the droplets do not wet the boundary, 
and that is why the semiperiodic systems behave essentially in the same way as periodic 
ones. For if s = 0.5, we expect the contact angle at this temperature to be about 0.7-7T 
and, therefore, boundary nucleation should prevail in small systems, as is also observed. 
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However, the crossover system size L x predicted in Eq. (51b) for these particular conditions 
is somewhat too small. We attribute this to the fact that the formula is expected to be valid 
for large systems, while our simulations are carried out for rather small L. Observations of 
configuration snapshots during our simulations reveal that the crossover to bulk nucleation 
takes place at larger L, between L = 20 and L = 100, without exhibiting a sharp crossover. 
In Fig. 8 we show snapshots of typical configurations with supercritical droplets in systems 
with L = 10 and L = 50. 

The switching behavior is rather different with zero applied surface fields, = 0.0. 
In that case, the droplets have a contact angle equal to n/2, which corresponds to a free 
boundary. Consequently, nucleation at the boundary is greatly preferred over nucleation in 
the bulk. As L increases, the system crosses over from the single-droplet regime into the 
multidroplet regime, but with the droplets still nucleating at the boundary. The maximum 
of the switching field in the vicinity of the thermodynamic spinodal is still visible, but the 
peak is much less pronounced than in cases with strong surface fields. This is due to two 
effects. First, the short waiting time makes the SD regime rather narrow, so that entropic 
terms in the free-energy barriers do not have much chance to assert themselves. Second, 
the entropic terms for nucleation at the boundary are smaller than those corresponding to 
nucleation in the bulk [see Eqs. (49) and (50)] and so is their effect on the switching field. 
This observation is in agreement with Eq. (53), which indicates that the crossover system 
size L x should diverge for 9 Y = n/2. For very large systems, the switching field increases 
and approaches the one measured for periodic systems. This is because the growth time 
becomes comparable to the nucleation time. The decay process then proceeds as follows. 
First, many droplets nucleate at the open boundaries and coalesce to create slabs of stable 
phase. Next, these slabs grow into the bulk until the magnetization is reversed. In this 
regime, to keep the lifetime constant while increasing L, one has to increase the field in 
order to compensate for the longer distance the growing interfaces have to travel. 

Finally, in Fig. 7 we see that Monte Carlo evaluations of the switching fields are in 
agreement with the theory discussed above. This makes these curves appear very similar 
to each other, even though the waiting times differ by a factor of thirty. Thus, we expect 
the theory to hold for all large values of r and to produce visually similar curves. This is 
important because a comparison of typical atomic time scales with typical experimental time 
scales leads us to believe that values of r relevant to experiments must be large numbers 
- much larger than can be simulated using Monte Carlo. In the next subsection, we will 
make theoretical extrapolations to this regime of large r for circular systems, which are more 
physically relevant than semiperiodic systems. 

Comparing the shapes of the switching-field plots in Fig. 7 with the theoretical predic- 
tions presented in Fig. 6, one can see that nucleation theory reflects the observed switching 
dynamics very well. 

B. Circular Systems 

Next we proceed with circular systems. Figs. 9(a) and (b) show the switching fields 
measured under the same conditions as those for the semiperiodic ones discussed above. 
One can see from the plots, as well as from simulation snapshots, that the behavior of 
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the circular systems is rather similar to that of semiperiodic lattices. There are two main 
differences. 

The first difference from the semiperiodic systems is that the switching field peaks are 
less pronounced. This effect is easy to understand when we compare the behavior of the free- 
energy barriers without entropy terms in Fig. 2. Since there is no sharp crossover between 
a coexistence regime and a single-droplet regime in circular systems, the maxima in the 
switching fields are shifted to larger system sizes and become flatter. In a system with weak 
boundary fields, like the one shown in Fig. 9(a), this can lead to the complete disappearance 
of the maximum. Instead, one observes a plateau which is eventually followed by an increase 
of the switching field at the crossover to the multidroplet regime. Although this is a typical 
situation, one can observe a weak maximum of the switching field even in circular systems 
with free boundaries, even without applied surface fields. This is demonstrated in Fig. 9(b) 
for i?s = 0.0. Though hardly discernible on the scale of the figure, there is a clear maximum 
in the switching field. The whole data set shown belongs to the single-droplet regime, and 
the decrease of the switching field would continue if we were to include even larger systems. 
Thus, the existence of the maximum of the switching field depends on the specific conditions 
under which the measurement is carried out. 

Second, one can observe small oscillations in the switching fields as functions of L for 
small systems. This is because in small systems the lattice structure in the vicinity of 
the boundary can influence the free-energy barriers (these types of corrections go beyond 
the theory presented in the previous section), and the interplay between the bulk and the 
boundary may depend significantly on the size of the system if it is small enough. 

C. Octagonal Systems 

The theory presented in Section IV treats the effects of the boundaries in a phenomeno- 
logical way, mimicking them by an effective interface tension assigned to that surface of 
the droplet which is also part of the system edge. Naturally, the structure of the lattice 
and interactions near the boundary can also influence the dynamics of the droplet surface 
fluctuations. These effects are not captured by our simple theory and may be expected to 
produce pre-exponential corrections to the metastable lifetime. To see whether our approach 
can really describe the essence of the physics in systems with different boundary structures, 
we have simulated the magnetization switching in octagonal systems with their boundaries 
modified as described in Subsection II A. 

For small systems on a square lattice, the best approximation to a circle is the same as 
the best approximation to an octagon. In order to observe a difference in shape, we need 
L > 20. For T = 0.9T C and r = 2000 MCSS, the switching field is nonzero only for L > 20 
(see Fig. 10). Even so, Fig. 10 shows that there is little difference between the switching 
fields for octagonal systems (shown in Fig. 1) and for circular systems (which have no special 
provision to ensure that boundary sites all have the same number of neighbors). 

Naturally, the difference between the octagonal and circular lattices becomes more dis- 
cernible at lower temperatures where the interfacial tension becomes anisotropic, and the 
role of the lattice irregularities at the boundary is also more important since the flipping 
probabilities of spins with different neighborhoods are further separated. However, our ob- 
servations indicate two things: 
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• For the weak fields at which a peak in H sw vs. L might be expected, the theory 
developed in the previous section for circular systems appears to be applicable to 
octagonal systems as well for the range of fields and waiting times employed in this 
study. 

• Geometric features much smaller than the critical droplet radius, such as omitted 
boundary bonds, are much less important than geometric features on the scale of the 
critical droplet radius, such as corners. 

Figure 11 shows how boundary fields change the qualitative dependence of the switching 
field H sw on the system size L of an octagonal system. For if s = 0, which corresponds 
to open boundary conditions without any modifications (i.e., Js = 0), we did not observe 
a peak in the switching field for these specific temperatures and waiting times. At both 
temperatures shown, the switching field monotonically increases with the system size and 
approaches its asymptotic value for very large systems. The absence of the peak is easy to 
understand in terms of the theory developed in the previous section: The crossover to the 
multidroplet regime, in which the switching field should be fairly constant, sets in before 
the decrease due to entropic terms can become discernible. As Hz is increased, the peak 
structure gradually develops and the switching field behavior becomes similar to that of a 
periodic system. One can see that for large enough systems the switching field is the same as 
that for periodic systems. The system size at which this occurs corresponds to the crossover 
from boundary nucleation to bulk nucleation as described in Subsection IV C. This behavior 
is clearly seen for a strong boundary field, H s = 1.0, at both temperatures. At T = 0.8T C , 
the droplets then practically do not wet the boundary, and even very small systems behave 
very much like periodic systems. For T = 1.3J the contact angle corresponding to H^ = 1.0 
is less than 7r and larger than n/2. In this case Eq. (26) is not expected to provide a good 
approximation for the boundary tension <r 2 and for the resulting contact angle 9y, because 
the boundary of the octagon with modified coordinations does not have the topology of a 
chain. Compared to circular systems, octagons are slightly less sensitive to the surface fields, 
which may be caused by the fact that not all their boundary spins interact directly. 

Similar effects, namely the disappearance of the switching-field peak, can also be observed 
for lattices which have modified spin-spin interactions close to their boundaries. This is 
shown in Fig. 12. In the switching- field behavior, these systems are more similar to the 
samples with H^ = than to samples with nonzero boundary fields. In addition to increasing 
the switching field, some irregularity is observed in the plot of H sw vs. L. This is apparently 
due to the fact that for small systems, the ratio of the length of a (10) face to the length of 
a (11) face depends nonmonotonically on L. This is a consequence of the fact that an ideal 
octagon cannot be constructed on a square lattice. The irregularity also indicates that there 
is still a difference between the (10) and (11) faces, in spite of the fact that boundary spins 
on both faces have two nearest neighbors. 

For large systems, Js > can be expected to yield similar results to H s > 0. This can 
most easily be seen for large systems with J 2 — > oo. With a single-spin-flip dynamic, it is 
impossible for either surface spins or the outermost layer of bulk spins to flip (see Fig. 1) 
in the limit oo. The coupling between these frozen spins and the unfrozen bulk spins 

makes the system very similar to an octagon of size L — 4 with H% = 1 on the (10) faces 
and with H^ = 2 on the (11) faces. Since the frozen spins still count towards the value 
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of m [Eq. (4)], the switching field will in this case be greatly increased for small systems. 
However, for larger systems, the behavior of the switching field is quite similar to other types 
of systems studied here. This is an expression of the fact that it is only the affinity of the 
boundary to the stable phase which is important for the switching behavior rather than the 
details of interactions in the vicinity of the system boundary. 

Thus, in spite of some quantitative differences between octagonal and circular systems, 
the basic behavior is very much the same, and we can conclude that the main features 
are well described by the nucleation theory from the previous section. On the other hand, 
it is seen that to achieve an accurate quantitative agreement between the theory and the 
simulation results, the theory would be required to include the dependence of the nucleation 
dynamics in the vicinity of lattice boundaries with different structures. 

VI. CONCLUSIONS 

We have studied the magnetization reversal dynamics in kinetic Ising systems with vari- 
ous kinds of boundaries. Our study emphasizes the stochastic regime of the space of magnetic 
fields and system sizes, in which the magnetization reversal is triggered by a single critical 
fluctuation. As expected, even minor modifications of the boundary conditions can lead to 
pronounced changes in the measured values of the metastable lifetime and switching fields. 
However, our simulations show that in spite of the quantitative differences, all the systems 
studied in this work are essentially described by the droplet nucleation theory outlined in 
Section IV. This theory allows us to understand in detail the behavior of the switching fields 
as functions of the temperature, system size, and waiting time. In particular, it can predict 
under which conditions it is possible to observe a maximum (peak) in the switching field as 
a function of L, and the extent of the single-droplet nucleation regime in the presence of 
boundaries. 

Magnetization switching in systems with boundaries displays a rich variety of nucleation 
and growth behaviors. Roughly speaking, the single-droplet and multidroplet regimes known 
from periodic systems come in two flavors in systems with boundaries. In each regime, the 
dominant nucleation sites can be located either close to the boundary or in the bulk of the 
system. Depending on the waiting time and the wetting angle, one can observe crossovers 
between these regimes, although the crossover regions can be rather wide. 

Qualitatively, the switching field behavior in the stochastic regime of a system with a 
boundary is similar to that in a periodic system if the contact angle of the droplets is larger 
than a right angle. If the contact angle is about n/2, the maximum in the switching field is 
in general suppressed and sometimes disappears completely. 

The structure of the system boundary can have a pronounced effect on the lifetime. 
Nevertheless, if the critical size of a droplet is larger than the irregularities at the boundary, 
then the switching field does not depend strongly on these details. The main contribution 
to the lifetime depends only on the attractive or repulsive influence of the boundary towards 
large droplets. However, for an accurate quantitative description of the switching dynamics, 
a complete understanding of the fluctuation dynamics of droplets close to the boundary may 
be required. 

Finally, we stress again that in this model the existence of a maximum in the switching 
field as a function of the system size has nothing to do with multidomain particles. The 
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systems studied here are single domain, and the switching field peaks are purely dynamic 
phenomena due to entropic effects and to a crossover between a regime in which the free 
energy of the critical fluctuation depends primarily on the applied field, and one in which 
it is determined by the system size. In systems with periodic boundary conditions this 
crossover is quite sharp and associated with a transition between well-defined single-droplet 
and coexistence regimes. Systems without periodic boundary conditions do not have a 
well-defined coexistence regime. However, the crossover behavior remains easily observable. 

Although the kinetic Ising model is certainly too simplified for quantitative agreement 
with experimental observations to be expected, it is tempting to ask to what extent our 
theory is capable of yielding results in a physically sensible range. For this purpose, we 
consider a hypothetical ferromagnetic monolayer which would correspond to our simulations 
of circular systems. We take its lattice constant as 0.3 nm, and its Curie temperature as 
375K. This implies a model exchange-interaction constant of J = 1.42 x 10~ 2 eV (compare 
with the exchange constant 1.19 x 1CT 2 eV of iron 111 ). With each atom we associate a 
magnetic moment of one Bohr magneton, /i B = 5.788 x 10~ 5 eV T _1 . This would correspond 
to S — 1/2 and a gyromagnetic ratio g — 2. Another necessary input is the relation between 
Monte Carlo time and real time, which is determined mainly by phonon frequencies. Thus, 
10 12 MCSS could roughly correspond to a physical time interval on the order of seconds. The 
prefactor which connects the exponential of the free-energy barrier in nucleation theory with 
the metastable lifetime is assumed to be of order unity. The switching fields for a waiting 
time of 10 13 MCSS, calculated for the above parameters at room temperature, comparable 
to 0.8T C (at which temperature the interface tension of the Ising model is a = 0.7459 J), are 
shown in Fig. 13. 

The simple insertion of physical parameters into our theoretical formulae, discussed in 
the previous paragraph, predicts switching fields that depend on the system size in a way 
very similar to experimental observations. However, the predicted switching fields are about 
one order of magnitude stronger than typical switching fields for nanoscale ferromagnetic 
systems. Also, the calculated switching fields grow with decreasing temperature or, equiv- 
alently, the calculated waiting times for realistic switching fields are unphysically long at 
lower temperatures. Nevertheless, keeping in mind the extreme simplicity of the model, we 
find this comparison very promising. Allowing for finite anisotropy of the model and includ- 
ing effects of disorder will decrease the switching fields significantly. Such investigations are 
already in progress. 9 The theoretical and numerical results of this study raise the exciting 
possibility that the convergence of more realistic microscopic spin models and further ex- 
perimental studies of the statistical switching behavior of well-characterized single-domain 
samples soon may provide quantitative, as well as qualitative, agreement between theoretical 
and experimental results. 
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TABLES 



TABLE I. The length I used in Eq. (16) to calculate -ffThSp- The listed value of r indicates the 
waiting time for which the Thermodynamic Spinodal was fitted to the peak in H sw vs. L. 



T/J 


r [MCSS] 


£ 


0.65 


10 7 


1.0(1) 


1.30 


4 x 10 3 


1.4(2) 


1.81535 


10 2 


2.0(4) 
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FIG. 1. An octagon of diameter L = 12. Surface sites, which contribute to the sum J2i s m 
Eq. (5), are circled. Surface bonds, which contribute to the sum Yl/ij)^ m Eq. (5), are shown as 
solid lines. In Eq. (3), all sites contribute to the sum J2i an d both boundary bonds and bulk bonds 
(dotted lines) contribute to the sum 
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FIG. 2. The normalized free-energy barriers (neglecting entropic corrections for the size de- 
pendent number of nucleation sites) for semiperiodic [AF n /AF§ D , solid curve, Eq. (37)] and cir- 
cular [AF /AFgj}, dotted and dashed curves, Eq. (48)] systems. The curve for periodic systems 
[AF I AFsd vs. |i7|/i?ThSp) Eq. (19)] coincides with the result for semiperiodic systems shown here. 
Note that the scaling variable \H\/ H^ hS differs from the scaling variable x = L/2R C in Eq. (48) 
only by a factor dependent on the contact angle 6y. Thus, with 6y and the magnetic field kept 
constant, these plots can be regarded as functions of the system size in units of the critical droplet 
diameter [see Eqs. (20), (38) and (48b)]. 
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FIG. 3. The contact angle By as a function of the boundary field Hy, for two different values 
of H and T. These plots are approximations calculated from Eqs. (26) and (27). Note that these 
(negative) values of H do not significantly influence the contact angle. This approximation gives 
an idea of the contact angles produced by the surface field in our Monte Carlo simulations in 
semiperiodic systems. Because of the different topology, this approximation may not be equally 
good for octagonal lattices. 
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FIG. 4. Three situations which may lead to the nucleation of a critical fluctuation in a semiperi- 
odic system: (a) a droplet at the boundary (the contact angle 9y [given by Eq. (27)] is shown), (b) 
the slab configuration, (c) a droplet in the bulk. The dashed lines represent edges of the system 
that are connected by periodic boundary conditions. The linear system size is L in both directions. 
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FIG. 5. The contact angle 6\ and the "angular diameter" 2(p of a droplet nucleating at the 
boundary of a circular system. R denotes the radius of the inner surface of the droplet, and L is 
the diameter of the circular system. 
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FIG. 6. An illustration of the theoretical results for the influence of 9y on the L dependence of 
the switching field. The solid and dashed curves correspond to surface nucleation; the dotted curve 
corresponds to bulk nucleation. For 9y = vr, bulk nucleation is dominant everywhere in the SD 
regime. For 9y = 3ir/4, surface nucleation dominates in the SD regime for L<L X , but for L>L X 
bulk nucleation becomes dominant. For 9y = n/2, surface nucleation dominates throughout the SD 
region. Note the following: (i) larger values of 9y produce more dramatic increases of H sw with L 
in the CE region; (ii) larger values of 9y produce more clearly pronounced peaks of H sw vs. L with 
larger values of H sw at the peak; (iii) the effect of bulk nucleation is to make the peak appear more 
pronounced. The curves were obtained by numerically solving for H using Eqs. (36), (49) and (50) 
for a fixed waiting time; both the parameters and the curves are in quantitative agreement with 
Fig. 7a. 
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FIG. 7. The switching field H sw as a function of system size L for various values of Hy, in 
semiperiodic systems at T= 1.3 J «0.57T C . Data for periodic systems are shown as dotted lines, 
(a) r = 1000 MCSS. (b) r = 30000 MCSS. Note the similarity in shape of the measured switching 
fields to the droplet-theory predictions presented in Fig. 6. 
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FIG. 8. Snapshots of Monte Carlo simulations for semiper iodic systems. Panel (a) shows a 
supercritical droplet nucleated at the boundary of a system with L = 10, and panel (b) shows a 
droplet nucleated in the interior of a system with L = 50. In both cases, T = 1.3 J, H% = 0.5 J, and 
for each system size, the magnetic field H is chosen equal to the switching field for r = 1000. Thus, 
the snapshots correspond to data points shown in Fig. 7(a). The smaller system (L = 10) is clearly 
in the "boundary dominated" single-droplet regime. The larger system (L = 50) is already in 
the crossover regime to "bulk dominated" nucleation, where the probability of observing a droplet 
nucleating at the boundary is smaller. For even larger systems with L « 100, almost all droplets 
nucleate in the bulk. 
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FIG. 9. The switching field H sw as a function of system size L for various values of Hy, for 
circular systems at T = 1.3J « 0.57T C . Data for periodic systems are shown as dotted lines, 
(a) r = 1000 MCSS. (b) r = 30000 MCSS. 
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FIG. 10. The switching field as a function of system size for small octagonal and circular 
systems with T = 0.9T C , r = 2000 MCSS, and H^ = J s = 0. 
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FIG. 11. The switching field as a function of system size for octagonal systems with various 
values of Data for periodic systems (dotted lines) were taken from Ref. 4. (a) T = 1.3 J w 0.57T C 
and r = 1000 MCSS. (b) T = 0.8T C wl.81J and r = 100 MCSS. 
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FIG. 12. The switching field as a function of system size for various values of Js for octagonal 
systems, (a) r = 100 MCSS, T = 0.8T C » 1.81J. (b) r = 1000 MCSS, T = 0.8T C » 1.81J. (c) r = 1000 
MCSS, T=1.3J «0.57T C . 



47 




FIG. 13. Droplet-theory predictions for the switching fields vs. system size for a hypothetical 
ferromagnetic monolayer with T c = 375 K, magnetic moment of a single atom of 1 Bohr magneton, 
and lattice constant 0.3 nm. According to Eqs. (3-5) in Ref. 7, such a system is expected to be 
single-domain if it is smaller than 400 nm. The three curves shown correspond to system boundaries 
with different properties characterized by the contanct angle 6y of the stable-phase droplets at the 
system boundary. We note that the crossover between the regime with droplets nucleating at the 
system boundary and the regime with bulk nucleation, as it was described in this work, is observed 
for the two upper curves. 
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